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ABSTRACT 

We present a detailed analysis from new multi-wavelength observations of the exceptional galaxy cluster 
ACT-CL J0102 4915, likely the most massive, hottest, most X-ray luminous and brightest Sunyaev-Zel’dovich 
(SZ) effect cluster known at redshifts greater than 0.6. The Atacama Cosmology Telescope (ACT) collaboration 
discovered ACT-CL J0102-4915 as the most significant Sunyaev-Zeldovich (SZ) decrement in a sky survey 
area of 755 square degrees. Our VLT/FORS2 spectra of 89 member galaxies yield a cluster redshift, z = 

0.870, and velocity dispersion, cr ga i = 1321 ± 106 km s _1 . Our Chandra observations reveal a hot and X-ray 
luminous system with an integrated temperature of T x = 14.5 ±0.1 keV and 0.5-2.0 keV band luminosity of 
L x = (2.19 ±0.11) x 10 45 /lygcrg s _1 . We obtain several statistically consistent cluster mass estimates; using 
empirical mass scaling relations with velocity dispersion, X-ray Y x , and integrated SZ distortion, we estimate 
a cluster mass of M200 = (2.16 ±0.32) x 10 15 We constrain the stellar content of the cluster to be less 

than 1% of the total mass, using Spitzer IRAC and optical imaging. The Chandra and VLT/FORS2 optical data 
also reveal that ACT-CL J0102-4915 is undergoing a major merger between components with a mass ratio of 
approximately 2 to 1. The X-ray data show significant temperature variations from a low of 6.6 ±0.7 keV at 
the merging low-entropy, high-metallicity, cool core to a high of 22 ± 6 keV. We also see a wake in the X-ray 
surface brightness and deprojected gas density caused by the passage of one cluster through the other from 
which we estimate a merger speed of around 1300 km s _1 for an assumed merger timescale of 1 Gyr. ACT- 
CL JO 102-4915 is possibly a high-redshift analog of the famous Bullet Cluster. Such a massive cluster at this 
redshift is rare, although consistent with the standard ACDM cosmology in the lower part of its allowed mass 
range. Massive, high-redshift mergers like ACT-CL J0102-4915 are unlikely to be reproduced in the current 
generation of numerical N-body cosmological simulations. 

Subject headings: cosmology: observations — galaxy clusters: general — galaxies: clusters: individual (ACT- 
CL JO 102^1915) — cosmic background radiation 


1. INTRODUCTION 

There are currently only a few examples of merging cluster 
systems in which there are spatial offsets (on the order of 200- 
300 kpc) between the peaks of the total and baryonic matter 
distributions. Some of these, like 1E0657-56 (the original 
“bullet” cluster at z = 0.296, Markevitch et al. 2002), Abell 
2146 (at z = 0.234, [Russell et al. l 2010), and possibly Abell 
2744 (at z = 0.308, |Merten et al.||2011) , contain in addition 
a cold, dense “bullet” of low entropy gas that is clearly as- 

t Based on observations made with ESO Telescopes at the Paranal Ob¬ 
servatory under programme ID 086.A-0425. 


sociated with the merger event. The offsets are due to the 
differing physical processes that act on the gas, galaxies, and 
dark matter. The gas behaves as a fluid, experiencing shocks, 
viscosity, and ram pressure, while the galaxies and dark mat¬ 
ter (or so we posit) are collisionless. These bullet systems 
have been used to offer direct evidence for the existence of 
dark matter ( |Clowe et al.||2004| |2006| ) and to set constraints 
on the self-interaction cross-section of the dark matter parti¬ 
cle ([Markevitch et al.|2004j [Randall et aL|2QQ8] [ Bradacet al.| 
|2008| ). Additionally, th e large merger velocity of 1E0657-56 
of around 3 000 km s _1 ( |Mastropietro & Burkert|2008[|Marke-| 
|vitch|2006| ), required to explain the morphology and temper- 
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ature of the gas from Chandra observations, is much higher 
than expected in cosmological simulations. These merging 
clusters provide interesting laboratories for several important 
topics in astrophysics and cosmology. 

Merging cluster systems can display a wide range of mass 
ratios, impact parameters, types of merging systems, and 
times since closest approach, and clusters evolve with cos¬ 
mic time. So it is worthwhile to find new bullet cluster sys¬ 
tems, especially at high redshift, to exploit their different ini¬ 
tial and final conditions and learn more about dark matter and 
the assembly of clusters through mergers. Massive clusters 
(~10 15 M q ) become increasingly rare at high z while thei r in¬ 
fall speeds are typ ically higher in the early Universe (see [Lee] 
& Komats u|2010| l. 

In this paper we r eport on ACT-CL J0102 -4915, a recently- 
discovered system ( |Marriage et al.|[20TT] ) which appears to 
be an excellent new example of a bullet cluster at redshift 
z = 0.87. It was found by the Atacama Cosmology Tele¬ 
scope (ACT, |Fowler et al.|2007|) through its strong Sunyaev- 
Zel’dovich signal (SZ, |Sunyaev& Zeldovich|1972|> , and was 
confirmed through optical and X-ray data by |Menanteau et| 
|aL[ ( |2010aT l. We demonstrate here that ACT-CL JO 102-4915 
is a rare and exceptional system at a time in cosmic history 
when the Universe was only half its current age. Addition¬ 
ally, it is a contender for being the most massive and X-ray 
luminous galaxy cluster at z > 0.6, likely competitor s be¬ 
ing CL J1226+3332 at z = 0.89 dMaughan , et al.|2004l> with 
M 200 = (1.38 ±0.20) x 10 15 /z^M© (fee & Tyson||2009]l and 
SPT-CL J21 06-5844 atz= 1. 14 with M 2 00 = (1.27 ±0.21) x 
1O 15 /i^M 0 ([Foley et al .[20 1 1|. Therefore, to reflect its excep¬ 
tional mass and to recognize the Chilean contribution to ACT 
we dub the cluster “El Gordo,” which means The Fat/Big One 
in Spanish. 

Throughout this paper we quote cluster masses as M 2 oo (or 
A/500) which corresponds to the mass enclosed within a radius 
where the overdensity is 200 (500) times the average (crit¬ 
ical) matter density, and we assume a standard flat ACDM 
cosmology with fl m = 0.27 and Oa = 0.73, and give relevant 
quantities in terms of the Hubble parameter H 0 = 1 Oh-10 km s -1 
Mpc -1 . For this cosmology the angular diameter distance to 
the cluster is 1617 Mpc and 1' corresponds to a physical 
scale of 470 hyp kpc at the cluster redshift. Magnitudes are 
reported in the AB system and uncertainties are reported for 
68% confidence intervals. 

2. OBSERVATIONS 

We have embarked on an ambitious multi-wavelength 
follow-up program aimed at characterizing the baryonic (i.e., 
gas, stellar) and dark matter mass components and the SZ- 
mass scaling relations of the most massive clusters discov¬ 
ered by ACT via the SZ effect. As part of this effort we 
have obtained total mass estimates based on galaxy dynam¬ 
ics (via measuring spectro scopic redshifts with 8-m class tele¬ 
scopes, |Sifon et al.|2011 [ 1, cluster X-ray emission with Chan¬ 
dra and stellar masses from infrared IRAC /Spitzer imaging. 
In the following sections we describe the different datasets 
and preview the results when applied to this study of ACT- 
CL J0102-4915. Images of the cluster from these various 
data sources are shown in Fig. [I] 

2.1. SZ observations 

ACT is a six-meter off-axis telescope designed for 
arcminute-scale millimeter-wave observations and it is situ¬ 


ated at an elevation of 5190 m on Cerro Toco in the Atacama 
desert in northern Chile (see |Fowler et al,|2010[ |Swetz et ah] 
|2011| for a complete description). One of the goals of this 
project is to measure the evolution of structure in the Uni¬ 
verse using massive galaxy clusters detected via the SZ effect 
(e.g |Sehgaletal.|201l| ). 

ACT observations began in late 2007 and to date it has 
surveyed two sky areas: one near declination -55 degrees 
(the southern strip) and the other on the celestial equator. 
The 2008 southern observations covered a region of 455 deg * 2 
and recently ACT has analyzed an additional 300 deg 2 along 
the equator (centered on the SDSS Stripe82, |Abazajian et al.| 
[2009l >. 

ACT-CL J0102-4915 was reported by ACT as a particu¬ 
larly strong SZ detection at a frequency of 148 GHz among 
23 high-significance c lusters from the 2008 southern survey 
(|Marriage et al.||2011). In Figure |T| (lower right panel) we 
show the filtered ACT 148 GHz intensity ma ps for ACT- 
CL J0102-4915 using a matched filtered (e.g., [Haehnelt &| 

|Tegmark|1996[[Melin, Bartlett, & Delabrouille|2006|l. 

The South Pole Telescope (SPT, |Carlstrom et al.|2009| l has 

reported its most significant cluster detections on a 2,500 
deg 2 survey ( [Williamson et al.||201l[ ), which overlaps with 
the ACT 2008 southern survey. This sample also contains 
ACT-CL J0102^1915; it is SPT’s most significant SZ clus¬ 
ter detection to date by nearly a factor of two, with a com¬ 
parable beam-averaged SZ decrement to the Bullet cluster 
1E0657-56. Among all clusters in the combined ACT and 
SPT sample, ACT-CL JO 1 02-4915 is exception al in having 
the highest SZ signal (ACT,|Marriage et al.|201 l]l or the high¬ 
est significance SZ detection (SPT, |Williamson et al.|2011) . 

2.2. Optical Imaging 

The initial optical observations of ACT-CL J0102-4915 
were carried out during December 9-12, 2009 using the SOI 
camera on the 4.1-m SOAR Telescope in Cerro Pachon (09B- 
0355, PI: Menanteau) using the SDSS griz filter set with ex¬ 
posure times of 540 s (6 x 90 s), 720 s (6 x 120 s), 2200 s 
(8 x 275 s) and 2200 s (8 x 275s) in g, r, i and z respectively. 
These exposure times were twice the nominal observations for 
the program as ACT-CL J0102-4915 was one of the targeted 
detections suspected to be a higher-redshift cluster. The run 
conditions were optimal and all four nights were photometric 
with seeing in the range 0"5 to 0"9. A complete de scription 
of the observation s and data analysis is provided in |Menan-| 
|teau et al7| ( ]2Q 1 Qa[ ). 

As part of our 2010B spectroscopic follow-up campaign 
of the southern cluster sample, we secured additional deep 
and wider Riz imaging (intended for slit mask design) and ex¬ 
tensive multi-object spectroscopy (MOS) with FORS2 on the 
Very Large Telescope (VLT, 086.A-0425 PLInfante). FORS2 
is a multi-mode (imaging, polarimetry, long slit and multi¬ 
object spectroscopy) optical instrument mounted on the UT1 
(Antu) telescope. We devoted a total of 12 hours of observ¬ 
ing time (service mode) on VLT to ACT-CL J0102-4915 of 
which 10 hours were devoted to MOS and 2 hours to imaging. 

In our original confirmation of ACT-CL J0102-4915 we re¬ 
ported the larger concentration of galaxies to the NW as the 
optical cluster counterpart to the ACT source and used the 
brightest galaxy it contained as the cluster center (see Fig¬ 
ure [T] top panels). However, even in the smaller » 4/5 x 4/5 
coverage of t he SOAR images (see Figure 9 from [Menanteau] 
|et al ■||20lTJat >, we noted a number of red galaxies with simi- 
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Figure 1. The multi-wavelength dataset for ACT-CL J0102-4915 with all panels showing the same sky region. (Upper left) The composite optical color image 
from the combined griz (SOAR/SOI) and Riz (VLT/FORS2) imaging with the overplotted Chandra X-ray surface brightness contours shown in white. The black 
and white inset image shows a remarkably strong lensing arc. (Upper right) The composite color image from the combination of the optical imaging from VLT 
and SOAR and IR from the Spitzer/IRAC 3.6pm and 4.5/rm imaging. The overplotted linearly-spaced contours in white correspond to the matched-filtered ACT 
148 GHz intensity maps. Bottom left False color image of the Chandra X-ray emission with the same set of eleven log-spaced contours between 2.71 cts/arcsec 2 
and 0.03 cts/arcsec 2 as in the panel above. The insert here shows the X-ray surface brightness in a cut across the “wake” region from the box region shown. 
Bottom right ACT 148 GHz intensity map with angular resolution of 114 and match-filtered with a nominal galaxy cluster profile, in units of effective temperature 
difference from the mean. The color scale ranges from -85 pK at the edges to -385 pK at the center of the SZ minimum. In all panels the horizontal bar shows 
the scale of the image, where north is up and east is left. 


lar photometric redshifts trailing towards the SE. The lack of 
coverage and large photometric redshift uncertainties made it 
difficult to confidently identify these galaxies as part of the 
cluster. 

The FORS2/VLT imaging aimed to obtain wider and deeper 
observations to confirm the substructure hint in the SOAR 
data and to secure galaxy positions for the spectroscopic 


observations. Our observations were carried out with the 
standard resolution of 0"25/pixel providing a field size of 
6/8 x 6/8 with exposure times of 640 s (8 x 80 s), 2200 s 
(10 x 220 s) and 2300 s (20 x 115 s) in R, I and z, respectively. 
The optical data were processed using a slightly modified ver¬ 
sion of the Python-based Rutgers Southern Cosmology image 
pipeline ( jMenanteau et al.|[2009| |2010b| l that takes us from 
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Figure 2. Histogram showing the redshift distribution for the 89 cluster 
member galaxies and the Gaussian with the same mean and width cr ga i. 
raw data to galaxy catalogs. Here we give a brief overview 
of the steps involved and refer the reader to the above pa¬ 
pers for a complete description. The initial standard CCD 
processing steps (i.e. flat fielding, bias subtraction) were per¬ 
formed via system calls to the ESO Recipe Execution tool 
(EsoRexrl In order to increase the depth and coverage for 
ACT-CL70102-4915 we co-added the imaging in the r, i and 
z and R, I and z from SOAR and FORS2 respectively into 
an iiber dataset. As the FORS2 RIz images were taken un¬ 
der clear but non-photometric conditions we utilize the SOAR 
photometric data to calibrate the photometry of the combined 
data. All science frames were aligned and median co mbined 
i nto a common pixel reference frame using SWarp ( |Bertin| 
|| |2006| l to a plate scale of 0"2525/pixel. Source detection 
and photometr y for the science catalogs were performed us¬ 
ing SExtractor ( |Bertin & Arnouts|1996| l in dual-image mode 
in which sources were identified on the iiber i-band images 
using a 1 ,5a detection threshold, while magnitudes were ex¬ 
tracted at matching locations from all other bands. 

2.3. Optical Spectroscopy and Redshifts 

Based on the possibility that ACT-CL J0102-4915 was ex¬ 
ceptional, we devoted the bulk (10 of 15 hrs) of the awarded 
time on VLT in semester 201 OB to securing spectra for galax¬ 
ies with colors consistent with being early-type c luster mem¬ 
bers. We refer the reader to a companion paper ( jSifon et ak] 
|2011] | where we fully describe the spectroscopic observing 
campaign of the 19 clusters at z > 0.3 in the ACT southern 
sample (of which ACT-CL J0102-4915 was part), including 
the data reduction, analysis and redshift determination. Here 
we summarize the most salient aspects of this work. 

The observations were executed in January 2011 under pho¬ 
tometric conditions with seeing of < 0"8. The observations 
cover the wavelength range ~ 4000-8000A using the GRIS 
300+11 Grism at 1.62 A/pixel and resolution A/AA ~ 660 at 
the central wavelength. At the cluster redshift of z = 0.87, the 
[O II] 3727A emission lines, the Call (K-H, 3950A) (which 
is the spectral signature of elliptical galaxies), plus other ab¬ 
sorption lines such as the G-band (4300A) as well as the 

2 http://www.eso.org/sci/software/cpl/esorex.htmt 


4000A break were located in the spectral range of the grism. 
These features facilitated the determination of galaxy red- 
shifts, galaxy cluster membership, and the galaxy velocity 
dispersion of the cluster. Integration times were 40 minutes 
per mask in order to maximize the number of galaxies while 
obtaining the necessary signal-to-noise ratio in the relevant 
lines. We observed a total of five FORS2/MXU masks in this 
configuration, providing redshifts for a total of 123 objects of 
which 89 were cluster members. 

The redshifts were measured by cross-correlating the spec¬ 
tra with galaxy spectral templa tes of the Sloan Digit al Sky 
Survey (SDSS) Data Release 7 (|Abazajian et al.|2009]» using 
the RVSAO/XCSAO package for IRAF (Kurtz & Mink|l998|). 
Most of these member galaxies are ellipticals with no emis- 
sion lines, and only a few emission line galaxies have been 
found to belong to the cluster. 

Cluster membership was determined using a rest-frame cut 
in velocity space of 40 00 km s -1 , and then ap plying the shift¬ 
ing gapper method of |Katgert et al.| ( |1996| l in order to re¬ 
move galaxies that lie outside a given gap in velocity space 
from the main body of galaxies. In particular, we defined 
the main body in bins of 15 galaxies and used a velocity gap 
of 500 km s -1 . We used the well-accepted biweight estimator 
( |Beers et al.|1990| l to estimate both the mean redshift of the 
sample, z = 0.87008±0.00010, and the velocity dispersiorj^of 
the cluster, cr ga i = 1321 ± 106 km s -1 . Errors were estimated 
via bootstrap resampling. The velocity dispersion of ACT- 
CL JO 1 02-4915 is larg er than all other clusters in the ACT 
sample (|Menanteau et al.||2010a| ISifon et al.||201 Iji and also 
larger than SPT-CL 2106 --5844 which is the largest in the full 
2,500 deg 2 SPT sample (Foley et al. 2011; Williamson et al. 
|2011| |. Using only the passive galaxies (i.e., objects with no 
emission lines) in ACT-CL J0102-4915, we obtained a mean 
redshift of z = 0.86963 ± 0.00012 with a velocity dispersion 
of 1280 ±108 km s” 1 . 

In Figure [2] we show the redshift distribution for all galaxy 
members from the FORS2 spectra. We see no indication of 
substructure along the line of sight, confirming the conclusion 
of |Sifon et al.| ( |2011| ) which used the |Dressler & Shectmaii| 
( |1988| ) test. However, as we will discuss in Sec. |3.2| we find 
substantial evidence for spatial substructure. 

2.4. Chandra X-ray Observations 

We were awarded 60 ks of Chandra observing time on this 
target during cycle 12 (PI: Hughes, Obsld: 12258) as part of 
our follow-up effort on the SZ clusters found in the ACT 2008 
455-deg 2 survey (|Menanteau et al.|2010a|. The Chandra ob¬ 
servation was carried out on January 26-27, 2011 for an ef¬ 
fective exposure of 59,260 s using the ACIS-I array. Figure [l] 
(bottom left panel) shows the surface brightness of the cluster 
after point source removal, exposure correction, and adaptive 
kernel smoothing. 

The cometary appearance of the X-rays is remarkable, ex¬ 
tending even to the apparent presence of two “tails” extend¬ 
ing off toward the NW. Orthogonal slices of the X-ray surface 
brightness show that these “tails” arise from significant de¬ 
pressions in the X-ray intensity, at a level of 20—40%, in a 
band extending from about 1' NW of the compact peak emis¬ 
sion region and continuing off to the NW out to the faint parts 
of the cluster. The depression is about 35" (270/t^' kpc) wide. 

3 As customary the quoted uncertainties for z and (T ga l only reflect statisti¬ 
cal errors from 89 objects and do not include systematic effects. 
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Figure 3. Chandra X-ray spectra of ACT-CL JO102-4915 for the total clus¬ 
ter in black, the cool bright peak (region 1 in Fig. [5) in blue, and the highest 
tempi l :gion (region 5) in red. The bottom panel shows the residuals 
between the best fit models (histograms in top panel) and data. 
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Figure 4. The X-ray bolometric l uminosity vs temp erature for a sample 
of well-studied clusters taken from |Markevitch| (T998) . The Bullet cluster 
(1E0 657-56) is the ope n square point at high temperature and luminosity 
from |Markevitchl(2006) . ACT-CL J0102-4915 is the red circle in the same 
region of the plot . 

We refer to this feature as the “wake”. 

The integrated Chandra spectrum is shown in Figure [3] 
(top black points). The extraction region for this spectrum 
excludes the central peak emission (from within the inner¬ 
most contour in Fig. [Tj and extends out to the outermost con¬ 
tour in Fig. [I] which is at a radius of around 2/2 (1 .0hj ( \ 
Mpc). An absorbed phabs*mekal model yielded a best- 
fit (source frame) temperature of kT = 14.5 ±1.0 keV, metal 



and the spectral extraction regions. (Lower panel) Best-fitted X-ray tempera¬ 
tures plotted versus spectral extraction region (see key in figure to left). Error 
bars are 1 -a. Region 1, corresponding to the bright peak, is the coolest region. 

abundance 0.19 ±0.09 with respect to solar values, and a 
bolometric luminosity of Lt, 0 i = 13.6 x 10 45 erg s -1 . Fig¬ 
ure!?] shows the L bo \-T x relation with ACT-CL J0102-4915 
added as the red point (note that a redshift co rrection fac- 
tor of (1+z) -15 was applied to L ho i following |Vikhlinin et| 
[al~| |2002[ >. The nearby point on this figure is the closest 
comparison cluster, 1E0657-56. Similarly, the X-ray lumi¬ 
nosity of ACT-CL J0102-4915 in the 0.5-2.0 keV band is 
Lx = 2.19 ±0.11 x 10 45 fc^erg s” 1 . 

The modest count rate of ACT-CL J0102-4915 (i.e., 0.2 s -1 
corresponding to a total number of detected X-ray photons of 
nearly 12,000) limits the number of regions into which the 
cluster can be subdivided for spatially resolved spectroscopy. 
We perform this measurement by dividing the total emission 
into roughly equal parts as shown in Figure [5] Spectra were 
extracted and fitted to the same absorbed thermal model as 
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Figure 6. The X-ray Chandra emission overlaid with the ACT 148 GHz data cont ours from Fig, IT] (left) a nd the predicted SZ signal at 148 GHz from the 
deprojected X-ray data (right). The black cross represents the SZ position reported by|Williamson et ar.|(2011) for ACT-CL J0102-4915. The contour levels are 
linearly spaced between the limits shown in brackets on each panel. The minimum central value ot the ACT and predicted signals agree to within ~15%. The 
white dotted line shows the symmetry axis chosen for the deconvolution. 


before; the best-fit temperatures are plotted vs. region num¬ 
ber in Fig. [5] The coolest (region 1, kT = 6.6 ± 0.7 keV) 
and the hottest (region 5, kT = 22keV) spectra are plot¬ 
ted in Figure |3] in blue and red, respectively. Note the strong 
Fe K line (observed at E = 3.56 keV) in the spectrum of 
the coolest region, which yielded a rather high metal abun¬ 
dance of 0.57 ± 0.20. Elsewhere the Fe K fine is not detected 
strongly, pointing to abundance variations in the cluster gas. 

We see that the peak emission region in ACT- 
CL J0102-4915 is an X-ray-bright, relatively cool, “bullet” 
of low entropy gas like in 1E0657-56. The enhanced Fe 
abundance suggests it was the core of a cool-core cluster 
that has, so far, survived a high velocity merger with another 
system. The cluster’s high overall gas temperature and spatial 
temperature variations (Fig. [3]) argue for significant amounts 
of shock heating of the gas. The steep fall-off in the X-ray 
surface brightness toward the SE, as well as the “wake” in 
the main cluster gas toward the NW, clearly indicate that the 
bullet is moving toward the SE. 

The two lower panels of Fig. p] show that the SZ and X- 
ray peaks are displaced by about T?2. Similar offsets between 
the peak SZ and X-ray emission have also been reported in 
the bullet-like cluster Abell 2146 at z = 0.23, for which Chan¬ 
dra and Arcminute Microkelvin Imager (AMI) observations 
at 13.9-18.2 GHz show a remarkably similar separation of 
the SZ and X-ray signal distributions (see|AMI Consortium:! 
|Rodrfguez-Gonzalvez et al.||2011| Figures 13 and 15). Two 
factors contribute to this. First, X-ray emission is more sensi¬ 
tive to denser parts of the cluster, since X-ray surface bright¬ 
ness is proportional to f n^T e a dl, where n e and T e are the elec¬ 
tron number density and temperature and a ~ 0.1 (for band- 
limited X-ray emission at high temperature), while the SZ sig¬ 
nal is proportional to the line-of-sight gas pressure / n e T e dl. 
Second, the combination of Chandra and ACT’s widely dif¬ 
ferent angular resolutions (< 1" vs. 1(4) is particularly rele¬ 
vant for observations of merging clusters with complex gas 
density substructure. 


In order to investigate these differences we developed a new 
technique for deprojecting X-ray images that, unlike tradi¬ 
tional methods utilizing azimuthally-averaged profiles, takes 
account of the full observed two-dimensional structure. The 
process assumes rotational symmetry about a user-defined 
axis and independently deprojects each column of pixels per¬ 
pendicular to that axis. For ACT-CL J0102-4915 we chose 
the symmetry axis so that it runs through both the peak X- 
ray emission at the bullet and the surface brightness depres¬ 
sion toward the NW. The axis is centered on R.A.=01:02:54.9, 
Dec.=-49:15:52.5 and oriented at a position angle of 136° 
(see Figure |6). Beyond a certain radius from the center 
(1/65 in the case of ACT-CL J0102-4915), where the sur¬ 
face brightness is background dominated, we used an analytic 
/3-profile fit to the azimuthally-averaged data to deproject the 
cluster emission smoothly to large radii. Interior to this ra¬ 
dius, using the adaptively-smoothed, 0.5-2.0 keV band im¬ 
age, we deproject the X-ray surface brightness in strips per¬ 
pendicular to the symmetry axis. Deprojection of the cluster 
emission is done in the usual manner (“peeling the onion”) 
using the best-fit integrated spectral parameters given above 
to relate the X-ray emission to the gas density of each depro¬ 
jected shell. We deproject the north and south parts of the 
cluster separately, but interpolate the density value across the 
cluster when “peeling” back any individual layer (to avoid a 
discontinuity at the symmetry axis). In general this process is 
not guaranteed to produce positive definite values (as physi¬ 
cally required), so we force the computed density to be zero 
in cases where the deprojection would result in a negative (or 
more precisely, imaginary) value. We quote electron density 
values using an electron to proton number ratio of n e /n p =1.2 
and for mass densities we assume a mean mass of 1 Am p . 

The deprojection results in an estimate of the three- 
dimensional distribution of the gas density in the cluster. 
From this we compute the predicted SZ signal at 148 GHz. 
Since the SZ depends sensitively on the gas temperature, as 
a first approximation, we use the fitted X-ray temperatures 
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in the 6 regions shown in Fig. [3] to generate a temperature- 
corrected SZ map of the cluster. This is then convolved 
with the ACT beam model ( jHincks et al.||2010) , added to 
a source-free region of the ACT map in the vicinity ACT- 
CL J0102-4915 (to include similar CMB and intrumental 
noise), and then is filtered with the same procedure used to 
detect ACT-CL J0102-4915. InFigure[6]we show contours of 
the predicted SZ signal overlaid on the Chandra data (right), 
next to an identical plot with contours of the ACT data (left). 
The centroid difference between the peak of the X-ray signal 
and simulated SZ map is reduced to 6", which is within the 
expected positional error circle. Furthermore the intensity of 
the predicted SZ signal agrees to within ~15% with the ACT 
data. To conclude, we note that the offset between the peak 
X-ray and SZ signals is predominantly a result of the differ¬ 
ent sizes of the beams between Chandra and ACT (at least a 
factor of 60). 


2.5. Spitzer/IRAC Imaging 

As part of a larger warm-phase Spitzer program (PID: 
70149, PLMenanteau) which aims to characterize the stellar 
mass content of clusters in the ACT sample (Hilton et al., in 
prep.), in August 2010 we obtained IRAC 3.6/irn and 4.5/jm 
observations of ACT-CL J0102-4915. Our observations were 
designed to provide coverage out to the virial radius for clus¬ 
ters at z > 0.4, using a 2 x 2 grid of IRAC pointings centered 
on the cluster position. A total of 10 x 100 s frames were 
obtained in each channel at each grid position, using a large 
scale cycling dither pattern. The basic calibrated data (BCD) 
images were corrected for pulldown using the s oftware of 
Ashby & Hora ^jand then mosaiced using MOPEX ( |Makovoz| 
|& Khan[2005l > to give maps which are « 13' on a side with 
0"6 pixel scale. The maps for each channel were then reg¬ 
istered to a common pixel coordinate system. By inserting 
synthetic point sources, we estimate the 80% completeness 
depths of the final maps to be « 22.6 mag (AB) in both chan¬ 
nels. 

Matched aperture photometry was performed on the IRAC 
maps using SExtractor in dual image mode, using the 3.6/jm 
channel as the detection band. We measure fluxes through 4" 
diameter circular apertures, which are corrected to estimates 
of total magnitude usin g aperture corrections (as measured 
by |Barmby et al.112008) of -0.35 ± 0.04, -0.37 ± 0.04 mag 
in the 3.6 pm, 4.5/im channels respectively. The photometric 
uncertainties were scaled upwards by factors of 2.8, 2.6 in the 
3.6pm, 4.5/xm channels respectively, in order to account for 
noise correlation between pixels introduced in the production 
of the mosaics which is not taken into account in the SEx¬ 
tractor error estimates. These scaling factors were d etermined 
using the method outlined in|Barmby et al.| (|2008ji. Finally, 
the uncertainties in the aperture corrections were added to the 
photometric errors in quadrature. 

3. ANALYSIS AND RESULTS 

3.1. Mass estimation 

Here we estimate the mass for ACT-CL J0102-4915 using 
all the mass proxies available from the X-ray data and optical 
spectroscopy. We present values in Tableland discuss them 
in detail in the next sections. In all cases we quote uncertain¬ 
ties that include measurement errors as well as the uncertain¬ 
ties in scaling relations between observables and cluster mass. 

4 See http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 

dataanalysistools/tools/contributed/irac/lixpulldown/ 


relation. For the dynamical and X-ray masses, scaling to the 
common mass estimator M 2 oo carries additional uncertainty 
that we quantify below and include in the Table. 


3.1.1. Dynamical Mass 

Here we report in more detail the dynamic al mass esti¬ 
mates for ACT-CL J0102-4915 calculated by |Sifon et ak| 
( |2011[ ). The velocity dispersion of the cluster is cr ga i = 1321± 
106 km s -1 . This matches the largest reported dispersion of all 
known clusters at z > 0.6 f rom RX J0152.7-13 57 at z = 0.836 
and (Tgai = 1322^g km s -1 (|Girardi et al.|2005|>, although note 
that the quoted velocity dispersion is tor a well-studied merg¬ 
ing system which is resolved into two components of con- 
sider ably lower ma s s. We use the M 200 -odm scaling relation 
from |Evrard et al.| (|2008|) to convert the measured velocity 
dispersion into a dynamical mass estimate, 


~ 0.7 h 70 (z) 


f <7DM \ 
\ ct DM,15/ 


l/« 

M 0 , 


( 1 ) 


where h-joiz) = h 2 oy/Cl m (l +z) 3 + Q ; \, ctdm.is = 1082.9 ± 
4.0 km s -1 , a = 0.3361 ± 0.0026 and <7 D m is the velocity dis¬ 
persion of the dark matter halo. The latter is related to the 
observed galaxy velocity dispersion by the velocity bias pa¬ 
rameter, b v = ggai/opM- The la test physically motivated simu¬ 
lations (see |Evrard et al.|2008l and references therein) indicate 
that galaxies are essentially unbiased tracers of the dark mat¬ 
ter potential, ( b v ) = 1.00 ± 0.05. Using a bias factor of b v = 1 
for the velocity dispersion for all galaxies, we obtain a dy- 
namical mass of M2oo,dyn = 1-86^0 49 x 10 15 fcjg Mr.,. Similarly, 
using only the passively evolving galaxies without emission 
lines, we obtain A^oo.dyn = 1.72^^ x 10 15 /f^oM 0 . We note 
that the |Evrard et al.j ( |2008| > scaling uses M 2 00 as the mass 
within a radius where the overdensity is 200 times the criti¬ 
cal densit y of the Universe, so the observed mass values re¬ 
ported in |Sifon et ah] ( |2011Jl have been scaled from critical 
to average density using = 1 - 16/rlm x This con¬ 
versio n factor was derived using a |Navarro, Frenk, & White] 
( |1997| ) mass profile (hereafter NFW) and the concentration- 
mass relation, c(M,z), from simulations (|Duffy et al.||2008Jl 
at z = 0.87 for the mass of the cluster. The reported uncer¬ 
tainties in the conversion factor reflect the <ri ogc = 0.15 scatter 
in the log-normal probability distribution of c(M,z). The ra¬ 
dius r2oo = 2111 ± 189 hjQ kpc was also calculated using M200 
and assuming spherical symmetry, using the average density 
of the Universe at redshift z, p(z) = il m (\ +z) 3 3// ( , / (87 tG). 


3.1.2. X-ray Mass 

We estimated the cluster mass fr om several X-ray mass 
pro xies using the scaling relations in[Vikhlinin et al.| 12009) 1 
and |Kravtsov, Vikhlinin, & Nagai| ( |2006| l. For the gas temper¬ 

ature and luminosity w e essentially follow the prescriptions in 
|Vikhlinin et al.| ( [20D9) > to apply those scaling laws to the X-ray 
data. We obtain masses of M 500 = 1 -39^q X 10 15 and 

M500 = 1.40^ x 10 1 5 hjQ Mq from the M 500 -T x and M 500 - 
L x relations, respectively. We refrain, however, from apply¬ 
ing the 17% "merging cluster mass boost" as implemented by 
|Vikhlinin et al.| ( |2U09) , since, given the exceptional nature of 
ACT-CL JO 102-4915, it is not obvious that such a correc¬ 
tion factor is appropriate. Due to the complex substructure of 
the cluster, for the M gas -based relations (Y x and M gas ), we use 
the 2D deprojection to generate an azimuthally-average radial 
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gas mass profile about the cluster center, which agrees well 
(within ~15%) with the azimuthally-averaged core-excised 
surface brightness profile. For the Y x mass proxy we obtain a 
mass of M 50 o = 1.55^15 x 1O 15 ^oM 0 . 

Since all of these mass proxies are potentially sensitive to 
merger boosts on T x and L x , we also explored using M gas , 
which is intrinsically unaffected by merger boosts. How¬ 
ever, the X-ray-derived measurement of M gas , which comes 
from the reconstructed density by deprojecting the surface 
brightness, is sensitive to clumping. This can in principle 
result in a bias, depending on whether clusters are indeed 
more highly clumped during mergers. In the case of ACT- 
CL J0102-4915, the excellent agreement between our pre¬ 
dicted SZ signal from the Chandra data and the actual mm- 
band ACT data suggests that gas clumping is not significant. 
We use the scaling law for Msotr-Mg as at redshift z = 0.6 from 
|Kravtsov, Vikhlinin, & NagaT| ( |2006| i, which yields a value of 
M500 = 1.67^20 x lO 13 ^ Mq and implies a gas mass frac¬ 
tion / gas = 0.133. The agreement among the four X-ray mass 
proxies is excellent: all are consistent within their l-cr ranges 
with a single mass value. 

Our M200 X-ray mass estimates are given in Table fusing a 
conversion factor of 1.85^ 2° from M500 to M2 qq. This factor 
is derived using the same procedure as in Section |3T.l| 

3.1.3. SZ Mass 

We estimate the SZ-derived mass of this cluster using the 
measured SZ flux of 490 ± 6O/2K as published in |Sehgal et| 
[ah] ( |2011| i. We combined this measurement with a fiducial 
scaling relation between y and M200 derived from simulations 
and calibrated by X-ray observations with best-fit parameters 
given in that same work. Using this scaling relation, we find 
the mass of this cluster to be ALjoo.sz = 1 x 10 15 H 2 qMq . 

This mass is consistent with the independently-derived mass 
reported by SPT ( [Williamson et al.|20~TT) i of M 2 00,sz = (1.78 ± 
0.64) x IO^/z^qMq, which uses a similar scaling relation be¬ 
tween the significance, andATjoo- This value also includes a 
6% reduction to the published mass value to reflect the change 
in redshift from their photometric estimate (z = 0.78) to the 
true spectroscopic value (z = 0.87). 

3.1.4. Stellar Mass from SED Fitting 

In order to measure stellar population parameters and the 
stellar mass content of the cluster using galaxy memb ers, 
we combine the optical photometry described i n Sec tion |2.2| 
with the IRAC photometry described in Section [23] by cross 
matching the catalogs using a 1"2 matching radius. We adopt 
the SExtractor MAG_AUTO measurements as estimates of 
total galaxy magnitudes in the optical bands, and the 4" diam¬ 
eter aperture corrected magnitudes as measurements of total 
flux in the IRAC bands. 

We fit the broadband Spectral Energy Distribut ions (SED) 
usi ng a si milar met hodology to that de scribed in |Shapley et| 
[ST] ( 120051 see also |Hilton et ~aI7||2010| >. We augmented the 
sample of 89 spectroscopic members to a total of 411 selected 
using photometric redshifts ( |Bemtez|2000| l re-calibrated with 
the 517 redshifts from seven clusters in the range 0.4 < z < 0.9 
from |Sifon et al.| < [201 1[ > and with the same griz photometry as 
ACT-CL .10102-4915. The improved photometric redshifts 
had typical Sz = 0.04 in the redshift range 0.8 < z < 0.9. We 
select galaxies with photometric redshifts within Az ± 0.06 of 
the cluster mean redshift._ 

We construct a grid of |Bruzual & Charlotj ( |2003[ > solar 



Figure 7. Examples of the SED fitting for four galaxy cluster members in 
ACT-CL J0102-4915. The solid line represents the best-fit SED for the 
broad-band g, r, i, z, 3.6/xm and 4.5/rm photometry shown as open symbols, 
metallicity models with exponentially declining star forma¬ 
tion histories with 20 values of r in the range 0.1 -20 Gyr, 
and 53 ages in the range 0.001-7.0 Gyr. We adopt a |Chabrier| 
( |2003| ) initial mass funct ion (IMF), and mode l the effect of 
dust extinction using the |Calzetti et al.| ( |2000| i law, allowing 
E(B-V ) to vary in the range 0.0-0.5 in steps of 0.02. The 
SEDs are fitted by analytically calculating the normalization 
N for each model SED where d\ 2 /dN = 0, adopting the model 
with the lowest x 2 value as the best fit. The stellar mass M* 
is then estimated from the value of N. We therefore fit for a 
total of four parameters (age, r, E(B-V), M*). Errors on the 
SED parameters were estimated using Monte-Carlo simula¬ 
tions. Figure [7] shows some examples of fitted SEDs. 

Given the small number of SED points compared to the 
number of free parameters, coupled with degeneracies be¬ 
tween several parameters (a ge, r, E(B-V); see the discus¬ 
sion in |Shapley et al.||2005] >, we do not expect the values of 
these parameters to be very robust. However, stellar mass is 
well constrained by the IRAC photometry, which probes the 
rest-frame near-IR at the cluster redshift, and stellar mass es¬ 
timates have been shown to be fairly insensitive to variations 
in other model p arameters, such as the assumed star forma- 
tion history (e.g.,|Forster Schreiber et al.|2004||Shapley et al.| 
|2005| . We note that uncertainties in the IMF and the modeling 
of thermally pulsating AGB stars (e.g., |Maraston|2005j |Con-| 
|roy et al.|2009| ) is likely to lead to the stellar mass estimates 
only being accurate to within a factor of two; these system¬ 
atic uncertainties are not taken into account in quoted errors 
on M*. The total stellar mass for all cluster members within 
r 2 00 =2111 kpc is M| 00 = 1.31 ± 0.26 x 1O 13 M 0 . This sug¬ 
gests a ratio of stellar mass to total mass of < 1% within 7-200 
for ACT-CL 10102^1915. 

3.2. Cluster Structure 

In the ACDM paradigm, the dark matter distribution should 
be aligned with the distribution of galaxies if the dark mat¬ 
ter particles are indeed collisionless. The cometary and 
bullet-like X-ray structure of ACT-CL J0102-4915 motivates 
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Table 1 

Dynamical Mass for ACT-CL J0102-4915 components 


Group IVgai 

redshift 

(km s 

(kpc hjo) 

Mi 00 

(1O 15 /^M 0 ) 

ACT-CL J0102-4915 (NW) 51 

0.86849 ±0.00020 

1290 ±134 

2062 ±233 

L76^;fi 

ACT-CL J0102-4915 (SW) 36 

0.87175 ± 0.00024 

1089 ±200 

1743 ±336 

L06^;g 

ACT-CL J0102-4915 (total) 89 

0.87008 ± 0.00010 

1321 ±106 

2111 ±189 

l.Stf&g 



Figure 8. The isopleths from the joint spectroscopically and photometrically selected galaxies in ACT-CL J0102-4915 using four methods. The isopleth contours 
using (upper left) the number density of galaxies per square arcminute, (upper right) the rest-frame optical t-band luminosity, (lower left) the rest-frame 3.6/j.m 
Spitzer/IRAC luminosity and (lower right) the SED-fitted stellar mass. The luminosity contours are both in units of 10 9 Lq arcsec -2 , and the stellar mass map is 
in units of 10 9 Mq arcsec -2 . The X-ray emission is shown as white contours using the same levels from Figure!!] All panels agree on the bimodal structure of 
the galaxy distribution; moreover, both the galaxy number and stellar mass density maps have peaks with statistically significant offsets with respect to the X-ray 
emission peak. 
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velocity (km/s) 



Redshift 


Figure 9. Redshift distributions for spectroscopic member galaxies in the 
NW (red) and SW (blue) clusters subcomponent. The underlying distribution 
for all of the 89 cluster member galaxies for ACT-CL J0102-4915 is shown 
in gray. 

our search for a spatial offset between the galaxies (that 
should primarily trace the dark matter distribution) and the 
X-ray gas in a cluster subcomponent as seen in, for example, 
IE 0657-558 dClowe et al.|2004), MS 1054-0321 dJeltema et| 
|al.|200ip , and MACS J 0025.4-1222 (jBradac et al.|2008| i. 

Our VLT campaign ( jSifon et al.|2011 [ i found no clear evi¬ 
dence for velocity substructure or a bimodal velocity distribu¬ 
tion in ACT-CL JO 102-4915. However, visual inspection of 
the optical images immediately suggested an elongated struc¬ 
ture in the distribution of the red galaxies along the NW-SE 
axis, coincident with the direction of the structure seen in the 
X-ray emission (see Figure [l] left panels). We investigate the 
spatial cluster structure in more detail by constructing four 
maps: a) the surface density of galaxies, b) the i-band lu¬ 
minosity density, c) the 3.6/xm/IRAC rest-frame luminosity 
density and d) the SED-fitted stellar mass density for spectro¬ 
scopically and photo metric ally selected cluster members as 
described in Section |3.1.4| These four maps use a 6"/pixel 
grid smoothed with a Gaussian kernel of width 20". Fig¬ 
ure [8] shows our deep optical multi-band imaging of ACT- 
CL J0102-4915 overplotted with yellow contours of the four 
density maps (one per figure panel) along with the same X-ray 
surface brightness contours from Figure [l] displayed in white. 

In all cases the galaxy distribution is obviously double 
peaked and a bright cD-type galaxy sits near each peak. For 
the galaxy density (Figure [8] top left) and stellar mass den¬ 
sity (lower right) maps, the peak of the X-ray emission lies 
between the peaks of the density. The luminosity densities 
for the i— and 3.6/zm-bands (Figure |8] upper right and lower 
left panels, respectively) are also double-peaked distributions. 
The i-band luminosity density, which probes the ultraviolet 
at z = 0.87, is somewhat unusual because it is dominated by 
the light contribution of the ultra-luminous Brightest Clus¬ 
ter Galaxy (BCG), discussed further below. As a result, the 
SE peak has the higher amplitude in the i-band, and coin¬ 
cides with the peak of the X-ray emission and with the lo¬ 
cation of the BCG. In the 3.6/zm band, which better traces 
the stellar mass, the NW peak is dominant, as with the galaxy 
number distribution. Note that the spatial distribution of the 
spectroscopically-confirmed cluster members is also double- 


peaked, similar to the contours in the top left panel of Fig¬ 
ure |8] 

Between and nearly equidistant from the two peaks in all 
four distributions lies a 20"-long straight blue arc, presum¬ 
ably the strongly lensed image of a high-redshift background 
galaxy (see black and white insert in Figure QY This feature 
appears in the trough between the two peaks of all four distri¬ 
butions. 

Regardless of whether we inspect the individual galaxies 
or their integrated light contribution, it is clear that there is a 
double peak structure in their spatial distribution with nearly 
parallel contours that also align with the direction of the lens- 
ing arc. We use this feature to define two components in 
the cluster (which we refer to as the NW and SE compo¬ 
nents) and investigate their masses and differences in veloci¬ 
ties of the components. Specifically, we divide cluster galax¬ 
ies into two groups along the line defined by the pair of RA 
and Dec positions (01:03:22.0, -49:12:32.9) and (01:02:35.1, 
-49:18:09.8). The velocity dispersions of the NW and SE 
components are 1290 ± 134 and 1089 ±200 km s _1 , respec¬ 
tively. We comp ute the ir masses using the same method de¬ 
scribed in section |3T±| and obtain 1.76^-® x 10 15 and 
1. 06^q 59 x 10 15 /z^qMq, respectively. Formally this implies a 
range of possible mass ratios from approximately 5:1 to 1:1 
between the cluster components; in subsequent analysis we 
will consider a 2:1 mass ratio. Table[l]gives computed values 
for the cluster components, as well as the cluster total. The 
sum of the masses of the two components determined sepa¬ 
rately, (2.8 ±0.9) x 10 i5 /?7qM q , is about l-cr higher than the 
total cluster mass inferred from the velocity dispersion of all 
galaxies. We discuss this more in the next section. 

The SE component has a line-of-sight peculiar velocity^Jof 
586 ± 96 km s -1 with respect to the main (NW) component. 
We also note that the cluster BCG (z = 0.87014 ± 0.00030), 
which is spatially located in the smaller SE mass component, 
has a peculiar velocity of 731 ± 66 km s -1 with respect to the 
main (NW) component. In Figure [9] we compare the distribu¬ 
tion of the velocities for the NW and SE components shown 
in red and blue respectively against the distribution for the full 
sample. 

We also use the cluster spatial information to obtain the stel¬ 
lar mass of the components wit h their respective ^oo value 
using the SED fits from section 13.1.41 Following the same 
component segregration as described above, we obtain stellar 
masses of (7.5 ± 1.4) x 1O 12 M 0 and (5.6±1.3) x 1O 12 M 0 for 
the NW and SW components respectively. 

3.3. Combined Mass 

Table [2] compiles our mass estimates for ACT- 
CL 10102^4915 from X-ray, SZ, and optical data sets. 
While central values range over nearly a factor of two, 
all estimates are statistically consistent. This can be seen 
by noting that the 2-<r mass ranges for all six mass esti¬ 
mates overlap for the mass range 1.78 x 1 0 i5 /!)q/W q and 
2.88 x 10 15 /j 7 qM q ; in other words, for any mass in this 
range, none of the six mass proxies are more than 2-cr away. 
For a combined mass estimate on which to base further 
conclusions, we construct a x 2 curve for a fit to a single 
combined mass value using the several mass estimates with 
their appropriate plus and minus error bars. For the three 

5 This is estimated as Av = c(zi -Z2)/(l +zi), the peculiar velocity of the 
SE component in the reference frame of the NW component. 
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Table 2 

Mass measurements for ACT-CL J0102-4915. 


Mass Proxy 

Measurement 

Scaling Law 

A/200 

(10 15 h~ 7 l o M 0 ) 

CTgal 

1321 ± 10610ns- 1 

M 200 -a 

1.86 4 q45 

Tx 

14.5± 1.0 keV 

A/500 -T x 

2.58+°J 4 

Lx 

2.19±0.11 

A/500 ~Lx 

2.6°! 1 ;|S 


3.20 ±0.24 

A/500 ~Yx 

2.88!°;™ 

A/gas 

2.2 ±0.1 x 10 14 M o 

A/500 -A/gas 

3-MSSg 

>’7 ’cmb 

490 ± 60 /rK 

A/200 -y 

l-64tl 

Combined 



2.16 ±0.32 


Note. — A/500 masses from Lx and Tx were scaled to M200 masses 
assumi ng an NFW density profile and the mass-concentration rela¬ 
tion of I Duffy et al~H2008) as described in the text. L x is in units 
of 10 45 /i^erg s _1 and is in the 0.5-2.0 keV band. Yx is in units of 
10 15 Mg keV. The combined mass is the minimum x 2 of the mass prox¬ 
ies in bold O ga i, Y x , and .yTcmb). 

independent mass estimators from velocity dispersion, Y x , 
and Sunyaev-Zeldovich decrement, we obtain a minimum x 2 
of 2.7 for 2 degrees of freedom (statistically acceptable) at a 
mass value of M 2 oo = (2.16±0.32) x 10 15 We do not 

combine all six mass proxies (which gives a somewhat higher 
mass value), since the four derived from the same X-ray data 
are all likely correlated to a degree that is difficult to evaluate. 

The above error range includes uncertainties from measure¬ 
ment error, mass scaling relations, and extrapolation to a com¬ 
mon M200 mass scale. Estimates of galaxy cluster masses 
also may be impacted by a number of potential systematic 
errors. Since cosmological conclusions are strongly sensitive 
to the estimated cluster mass, it is important to have a realistic 
appraisal of these errors. Here we describe several relevant 
sources of systematic error or biases in both measured mass 
estimates and in the theoretical comparison with cosmological 
models. 

Dynamical mass estimates are based on the assumption that 
the system is dynamically relaxed, which is almost surely not 
the case in ACT-CL J0102-4915, based on the observed spa¬ 
tial galaxy distribution. If we are observing a merger with the 
component velocities primarily in the plane of the sky, then 
significant kinetic energy can be contained in motions trans¬ 
verse to the line of sight which are not observed in galaxy 
redshifts; the dynamically inferred mass would be lower than 
the actual system mass. This seems to be consistent with the 
measurements of the dynamical masses of the two subcom¬ 
ponents in ACT-CL .10102-4915 (explored in Section |3.2[ >. 
which when summed, (2.8 ±0.9) x 10 15 is in better 
agreement with the X-ray-derived masses than the full galaxy 
sample. Therefore the dynamical mass of the full cluster 

quoted in Table [2] may be an underestimate. _ 

A number of groups (e.g.,|Ricker & Sarazin|2001||Randall| 
|et al.||2002)> have shown that cluster mergers can result in a 
transient boost in X-ray temperaure and luminosity during a 
merger, compared to the X-ray emission of the merged clus¬ 
ter once it relaxes. The effect that such boosts have on the 


X-ray mass proxies has been more r ecently investigate d by 
|Kravtsov, Vikhlinin, & Nagai|P006Jl,|Poole et al.|(j2007j) and 

|Yang et al.| ( |2010| >, among others. |Vikhlinin et al.| ( |2009b ar¬ 

gue that masses estimated from T x should be increased by a 
factor of 1.17 (which we do not apply) for merger systems, to 
account for the fraction of merger kinetic energy that re mains 
unthermalized. |Kravtsov, Vikhlinin, & Nagai| ( |2006| l intro¬ 
duced the Y x mass proxy which has been shown to be fairly 
independent of merger dynamics (e.g.,|Kravtsov, Vikhlinin, &| 
|Nagai|2006| |Poole et al.|2007] l. ACT-CL JO 102-4915 shows 
strong evidence for being in the merger process, particularly 
the compact region of X-ray enhancement and the double- 
peaked sky density distribution of galaxies. So in our com¬ 
bined mass estimate we strongly prefer to use Y x , but we note 
that the masses for all the X-ray proxies are consistent with 
each other. In particular the agreement of our masses from 
the Y x and M gas proxies suggest that our X-ray masses are not 
significantly boosted by the merger. 

The accuracy of the mass scaling law for the SZ proxy (here 
>’7cmb) is yet to be determined and in fact is one goal of the 
ACT follow-up program. Still it is interesting to note that the 
inferred SZ mass for ACT-CL J0102-4915 is only 60% of the 
Lc-inferred mass, albeit consistent within the large errors. The 
excellent agreement we find between the predicted SZ signal 
from the Chandra and ACT data (as illustrated in Figure [6ji 
suggests that we have a realistic gas model for the cluster and 
thus points to the mass scaling law, particularly at high red- 
shift, as the main source of possible discrepancy. A similar 
difference between the SZ and X-ray derived masses is also 
seen for the massive cluster SPT-CL J2106-5844 at z = 1.132 
( |Foley etal.|2011| ), while |Andersson et aT7| ( |2010[ ) find that SZ- 
derived masses are some 22% less than their X-ray masses. 

Given six mass proxies which are all statistically consistent, 
it seems unlikely that our combined mass estimate is subject 
to significant systematic error. The known systematics dis¬ 
cussed here tend to move central mass estimates towards even 
closer agreement - a possible modest increase in the mass es¬ 
timated using the galaxy velocity dispersion and the SZ mass 
scaling relation, and a possible modest decrease in the mass 
estimated from X-ray scaling relations, apart from Y x . 

For cosmological model comparisons, the theoretical mass 
function for clusters this massive is not precisely determined. 
In principle, direct numerical simulations of dark matter cos¬ 
mology in large volumes and low mass resolution can de¬ 
termine the high end of the mass function in a straightfor¬ 
ward way. The largest simulations to date have comoving 
volumes of roughly 2000 Gpc 3 , on the order of our Hub¬ 
ble volume, which nonetheless contain perhaps only a sin¬ 
gle cluster with mass 10 15 M Q at z = L An analytical model 
of the high-mass mass function accurate at the 10% level 
has not yet been demonstrated, so the mass functions used 
for comparison involve extrapolation from lower-mass halos 
where statistics from numerical simulations are better. Care 
also must be taken to ensure that initial conditions in simula¬ 
tions are sufficiently accurate; the commonly used Zeldovich 
Approximation for initial conditions may underestimate the 
number of high-mass haloes by 10% to 30% compared to 
more accurate second-order Lagrangian perturbation theory 
(|Crocce, Pueblas, & Scoccimarro||2006l |Berlind, Farrar, &| 

|Zaw|2010| i. Finally, large-volume simulations give dark mat¬ 

ter halo masses, which then must be translated into total halo 
masses, with possible systematic departures from the com¬ 
monly assumed universal baryon to dark matter ratio in large 
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Figure 10. The optical color-magnitude relations for ACT-CL .10102-4915 
for spectroscopically confirmed cluster members. Absorption and emission 
line galaxies are represented by red circle and blue squares respectively while 
the BCG is denoted by a yellow circle. The one spectroscopically confirmed 
AGN is represented by a cyan star. The dashed lines represent the fitted rela- 
tions for the cluster RXJ0152-1357 at z = 0.83 from|Blakeslee et al.|<2006) 
for the corresponding HST/ACS filters redshifted to z = 0.87 1'he gray tri¬ 
angle represents the observered magnitude of a M* galaxy (i = 22.16) in the 
CMR. 
clusters. 

3.4. Cluster Evolution 
3.4.1. Color-Magnitude relations 

Considerable effort has been devoted to understanding the 
formation and evolution of ellipticals in clusters up to z — 1.4 
(e.g.,|Lidman et al.|2008||Hilton et al.|2009). A consistent pic¬ 
ture has emerged in which these galaxies are a well-behaved 
class of objects with structural and chemical properties obey¬ 
ing simple power-law scaling relations. The high degree of 
homogeneity observed in these relations and in their stellar 
population suggests formation at hi gh redshift followed by 
passive evolution of their galaxies (Bower, Lucey, & Ellis 

119921 IBarrientos, Schade, & Lopez-Cruz||1996l lEllis et al. 

|l997[|Holden et al.|2005||Mei et al.|2009| l. In principle, when 

precise photometry is available measurements of the intrinsic 
scatter about the color-magnitude relation can be used to con¬ 
strain the ages of the constituent stellar populations of cluster 
galaxies. 

Here we study the optical and optical/IR color-magnitude 
relations (CMR) for ACT-CL J0102-4915 to assess the evolu¬ 
tionary stage of the cluster. We construct the color-magnitude 
diagrams for spectroscopically confirmed cluster members us- 


Figure 11. The optical to infrared color-magnitude diagram for cluster mem¬ 
ber galaxies. Symbols are as in Fig.Qo] The dashed lines aga in represent the 
fitted relations for the cluster RXJOlsTll 357 at z = 0.83 from IBlakeslee et al.l 
(2006) , now transformed to the Spitzerl IRAC bands as described in the text. 

ing the basic approach described in|Blakeslee et al.| (|2006J». 
For total magnitudes we use the S Ex tractor MAG_AUTO 
values (with a Kron factor of 2.5) and we measure colors 
within the galaxy effectiv e radii R e which w ere computed us¬ 
ing GALFIT (version 3.0, |Peng et al.|2010) to fit each galaxy 
to a |Sersic| ( fl9681 > model from the PSF-convolved z'-band im¬ 
ages. Bright nearby objects are fit simultaneously while faint 
ones are masked. Magnitudes in each passband are then mea¬ 
sured within a circular aperture of radius R e while enforcing 
a minimum radius of 1.5 pixel. 

In Figures [TO] and |TT] we show the optical and optical-to- 
infrared color-magnitude diagrams for spectroscopically con¬ 
firmed cluster members. Absorption and emission line galax¬ 
ies are represented by red circles and blue squares respectively 
while the BCG is denoted by a larger yellow circle. In order to 
ease the comparison between panels all colors are shown with 
respect to the total z-band magnitude. In order to compare 
with the expected colors of passively evolving system in fig- 
ure[l0]we show for comparison the fitted optical CMRs for the 
cluster RXJ0152-1357 at z = 0.83 from |Blakeslee et al.l l |2006) 
which was observed with the corresponding HST/ACS filter 
set and red s hifted to z = 0.87 using a solar metallicity |Bruzual| 
|& Chariot] (|2003)> burst model with r = 1.0 Gyr formed at 
Zf = 5. We also display for reference the typical magnitude 
for an M* (z = 22.16) galaxy at the cluster redshift in the CMR 
as a gray triangle. In Figure [Tl] we show the fitted i—z vs. i 
relation from Fig. [10] transformed to the SpitzerflRAC band 
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magnitude as a function of redshift for elliptical galaxies over a wide range 
intrinsic luminosities (i.e. the curves represent the observed magnitude for 
galaxies of the same luminosity). 

using the same stellar population model. 

Despite some scatter, a “red sequence” of passively evolved 
galaxies is quite distinguishable and it is in good agreement 
with the observed relations for RXJ0152-1357. However, un¬ 
like IBlakeslee et al.| (|2006|> we see bluer early-type galaxies 
brighter that M~ possibly due to recent star-formation trig¬ 
gered by the merger of the cluster. We note however, that our 
classification of galaxies is only based on the color and visual 
information in our seeing-limited ground-based imaging. 

From the figures is also clear that the BCG is extremely 
luminous (~ 1.5 mag brighter that any other galaxy) and sig¬ 
nificantly bluer in its optical colors than expected from a pas¬ 
sively evolving elliptical galaxy. In the following section we 
describe the BCG properties in more detail. 


3.4.2. Luminosity of Brightest Central Galaxy 

In Figure [12] we plot the observed z-band magnitude of 
the BCGs ofACT clusters in the southern sample, plus the 
BCG’s for three X-ray-selected, massive, high redshift clus¬ 
ters: RDCS 1252.9-2927 dBlakeslee et al.|2003||Rosati et al.| 
[2004l> at z = 1.237, XMMU J2235.3-2557 ([Mullis et al.|2005[ 
[Rosati et al. |2009) atz = 1.393 and XMMXCS J2215.9-1738 
( |Stanford et aLp 006; |Hilton et al.|[2009] > at z = 1-457. For 
clarity the BCG for ACT-CL 10102^1915 is shown as the yel¬ 
low star, with a luminosity of ~ 22 L* in the z-band. This is 
clearly an extremely luminous object and it is the brightest 
BCG in the ACT southern sample. We also plot the expected 
apparent z-band magnitudes of elliptical galaxies of differ¬ 
ent luminosities as a function of redshift. We use L* in the 
z'-band as defined by |Blanton et al.| ([2003]! at z = 0.1 (M* = 
-20.82 + 51og(0.7 )) and allow passive evolu tion according to 
a solar metallicity,|Bruzual & Charlot|d2003[! r = 1.0 Gyr burst 
model formed at z/ = 5. The resulting magnitude-redshift re¬ 
lation is plotted for a range of luminosities (0.4L*, L *, 4L*, 
8L* and 16 L*) aimed at representing the brightest observed 
BCGs. 

The VLT/FORS2 spectrum of the BCG in ACT- 


CL J0102-4915 has the signature of an E+A+[0 II] galaxy. 
This galaxy is a possible high-redshift analog of NGC 1275, 
the blue a nd radio-loud cD galaxy in the Pe rse us cl uster at 
z = 0.018 dMcNamara, O’Connell, & Sarazin:j|1996| l. Blue 
BCGs like this one usually live in the centers of cool core 
clusters, and in this case the BCG is located within the cool, 
low entropy, merging core in ACT-CL J0102-4915 revealed 
by the Chandra data. We also obtain the physical size of the 
BCG from the structu ral parameters estimated using GAL- 
FIT in section |3.4.1| The galaxy has an effective radius 
R e = 10.4 /zjfjkpc for a Sersic index n = 1.3. It is likely that 
this BCG is the result of a recent galaxy merger; we speculate 
that this extreme galaxy resulted from a merger of two BCGs 
during the cluster collision. 

4. DISCUSSION 
4.1. A New Bullet Cluster? 

Given the absence of a weak-lensing mass reconstruc¬ 
tion, we use the galaxy distribution as a proxy for the to¬ 
tal mass distribution. This should be a robust approxima¬ 
tion since for both 1E0657-56 dRandall et al.||2008|l and 
MACS J005.4-1222 ([Bradac et al.|2008(i the galaxy and to¬ 
tal matter distribution closely track each other. Figure |8] (up¬ 
per left panel) shows that the SE galaxy distribution precedes 
the peak gas location by ~ 22" (173 Ajgkpc) in the approxi¬ 
mate direction of motion of the merger. This suggests ACT- 
CL J0102-4915 as another “bullet” cluster undergoing a high- 
velocity merger event resulting in a significant spatial sep¬ 
aration in the main baryonic (gas) and total mass compo¬ 
nents. Moreover, ACT-CL JO 102-4915 is at a significantly 
earlier cosmic epoch, a lookback time around 4 Gyr longer 
than 1E0657-56. However, this picture requires confirmation 
through weak lensing mass mapping to determine the spatial 
distribution of the dark matter. 

4.2. Rarity of ACT-CL J0102-4915 

ACT-CL J0102-4915 is a rare, massive high-redshift clus¬ 
ter. We use the convenient fitt i ng for mulas provided in ap¬ 
pendix A of |Mortonson et al.| ( |201 1[ > to compute the con¬ 
fidence level at which the existence of a cluster like ACT- 
CL J0102-4915 would be expected in standard ACDM cos¬ 
mological models, where bound objects grow via gravita¬ 
tional instability from Gaussian initial density fluctuations. 
In Figure [13] we show AJioofe) exclusion curves for which a 
single cluster with mass M 2 oo above the corresponding curve 
would conflict with ACDM and quintessence at 95% confi¬ 
dence level, including both sample and cosmological param¬ 
eter variance. In other words, the exclusion curves represent 
the mass threshold as a function of redshift for which a cluster 
is less than 5% likely to be found in a survey region for 95% 
of the ACDM parameter variance. 

For purposes of addressing how rare ACT-CL J0102-4915 
is, in Fig. [13] we plot the exclusion curves for the full sky, the 
region analyzed for the ACT survey (755 deg 2 ) where ACT- 
CL J0102-4915 was discovered, and the combined ACT+SPT 
survey region (2800 deg 2 of which 455 deg 2 were observed 
by both ACT and SPT) as this is the most massive cluster at 
z > 0.5 in the combined area and detected by both experi¬ 
ments. Our combined estimated 1 -a mass range for the clus¬ 
ter (in Tabled) ranges from 1.84 to 2.48 x 10 15 hflM & . At the 
cluster redshift of 0.87, the lower end of this mass range falls 
below the ACT+SPT region exclusion curve, while the upper 
end of this range falls well below the full-sky exclusion curve. 
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parameters variance. The black dashed and dotted lines represent the exclu¬ 
sion curves for the full sky and the ACT survey analyzed region of 755 deg 2 
(i.e. the mass limit for a cluster at a given redshift that it is less than 5% 
likely to be found in 95% of the ACDM parameter probability distribution). 
We also show the same exclusion limit for 2800 deg 2 (ACT+SPT area) as the 
black solid curve. The open circle shows our combined mass estimate for 
ACT-CL J0102-4915, with l-a error bars. 

In other words, the cluster is not an unlikely occurrence in the 
ACT+SPT sky region, provided that its actual mass is l-cr or 
more below our nominal mass of 2.16 x IO^/z^L/q. If the 
cluster mass is 1 to 2 -a above this nominal mass, then it is 
a rare occurrence in the ACT+SPT sky region but not unex¬ 
pected in the entire sky. The cluster mass would have to be 
about 3-cr above the nominal mass for it to lie above the 95% 
exclusion curve for the entire sky. We conclude that, while 
this cluster is clearly rare, it is not massive or early enough to 
put significant pressure by itself on the standard cosmological 
model. 

Because the high-mass cluster mass function is steeply 
falling with mass, the signals we observe are more likely 
to be from a cluster with mass below our nominal mass of 
2.16 x 10 15 AyQ A/q and a positive measurement error than 
from a cluster with mass above this level and a negative mea¬ 
surement error. In other words, the cluster mass is more likely 
to be near t he lower end of our nominal mass range than the 
upper end. |Lima & Hu] ( |2005| l give a discussion of this Ed¬ 
dington bias and compute a correction for this effect, which in 
our case is of order 5%. This correction does not qualitatively 
c hange any conclusion s presen ted here. 

|Lee & Komatsu] ( |2010| l investigated the rarity of 
IE 0657-558 in cosmological simulations down to z < 0.5 
based on its large merger velocity and found the existence 
of such a merging cluster to be in tension with the expecta¬ 
tions of ACDM. They suggested that such high merger speeds 
could potentially be more common at high redshift, while, on 
the other hand, massive, ~ 1O 1 5 M 0 clusters are rare at ear- 
lier times. Using simulations, |Forero-Romero et al.| (|2010]> 
studied the expected probability distribution of the displace- 
ments between the dark matter and gas cores, as observed in 
IE 0657-558, and found it to be expected in l%-2% of the 
clusters with masses larger than 10 14 M Q . 


Here we briefly consider the rarity of ACT-CL J0102-4915 
as a massive cluster undergoing a major merger with a mass 
ratio of ~ 2 to 1 at z = 0.87, ignoring its merger velocity and 
spatial separation between gas and dark matter. We exam¬ 
ined the output of a large-volumne (3.072/r 1 Gpc) 3 ACDM 
cosmological N-body simulation with 29 billion particles per¬ 
formed by I. Iliev with the Cubep3m code at NIC Juelich 
(G. Yepes, private communication). Within this large vol¬ 
ume there is only one cluster of mass 1.9 x 1O l5 fcjgM 0 at 
z — 1 , which implies a very low number density of around 
10 -11 Mpc -3 for such clusters. This system has only one small 
potential substructural feature with a mass of around 10 13 M Q , 
suggesting a recent minor merger with mass ratio of 100 to 
1. The comoving volume of the entire sky between z = 0 and 
z = 1 is 54.7 hr 3 Gpc 3 and hence the examined simulation sam¬ 
ples roughly half of the total cosmic volume accessible by SZ 
cluster observations. We conclude that a merger event of the 
mass and mass ratio similar to the one we are witnessing in 
ACT-CL J0102-4915 is a rare event in the current set of large 
volume simulations. Significantly larger simulations, or im¬ 
proved analytic approximations to halo growth and evolution, 
will be required to assess directly the likelihood of massive 
high-redshift merger events such as ACT-CL J0102-4915. 

4.3. Merger Speed 

We have argued above that the depression in the X-ray sur¬ 
face brightness northwest of the merging bullet is a “wake” 
caused by the passage of one cluster through the other. In 
order to study this feature in more detail, we examine the 
map of the deprojected electron density in the cluster mid¬ 
plane (i.e., the plane of the sky) of ACT-CL J0102-4915 (see 
Figure [14]! . The merging bullet is indeed dense with a peak 
density of n e = 0.047 cm -3 , while the outskirts fall to a value 
n e ~2x 10 -4 cm -3 . Formally our simple deprojection proce¬ 
dure fails in the region of depressed surface brightness (pre¬ 
sumably because the cluster is not axially symmetric as as¬ 
sumed) and the density values here were pegged to zero; this 
area appears black in Fig. [14] 

Nonetheless, clearly the NW region of depressed X-ray sur¬ 
face brightness requires a region of extremely low density in 
the cluster interior. The region of depressed density is about 
0.8-1 ' (370/Z7Q -460/zyQ kpc) wide. It begins approximately 
2 . 6 ' (1200/z7o kpc) behind the leading edge of the bullet to¬ 
ward the northwest and can be traced into the core of the clus¬ 
ter to within about 1' (460^ kpc) from the bullet. At locations 
even closer to the bullet, the deprojected density map shows a 
complex, disorganized set of peaks and valleys suggesting the 
presence of merger-driven turbulence behind the cold merging 
core. 

For how long after the passage of a merging bullet should 
we expect a wake to persist? The natural timescale for the 
re-establishment of pressure equ ilibrium after a merger is the 
cluster sound crossing time (e.g., |Ricker & Sarazin|2001| ). As¬ 
suming a constant mean cluster gas temperature (about 14 
keV), the sound speed is c s ~ 1900 km s -1 , which implies 
a sound crossing time to the virial radius (here R200 —2.1 
Mpc) of t sc ~ 1 Gyr. If we assume that the displaced gas in 
the wake returns back to refill the cavity on this timescale, 
then we can make an independent estimate of the speed of the 
merger. Since the wake persists for a projected distance, d, of 
1200/zyo kpc behind the bullet, we estimate the merger speed 
to be v merg — d/t sc /cos(0) = 1200/cos (6) km s -1 , where 6 is 
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Figure 14. Map of the electron density in the mid-plane of ACT- 
CL J0102-4915 from a deprojection of the Chandra image. The contours are 
at density values of 0.002, 0.003, 0.004, 0.006, 0.008, 0.011, 0.016, 0.023, 
0.032,0.045 cm- 3 . 

the inclination angle of the merger direction (with respect to 
the plane of the sky). The merger speed can also be estimated 
from the peculiar velocity difference between the SE and NW 
galaxy concentrations as v mer g = v pec / sin(0) = 600/ sin(0) km 
s -1 . With these values the inclination angle becomes 6 = 21° 
and the merger speed is 1300 km s -1 . If the gas falls back 
into the cavity wake more quickly, that would decrease the in¬ 
clination angle and increase the calculated merger speed. If, 
by way of illustration, we halve the refilling time to 0.5 Gyr, 
then the inclination angle becomes 6 = 14° and the merger 
speed increases to a value of 2400 km s -1 . Clearly a more 
accurate estimate of the merger speed will require detailed 
N-body/hydro simulations, but a rough picture of the merger 
dynamics in ACT-CL J0102-4915 emerges from basic argu¬ 
ments about the wake morphology in the system. 

5. SUMMARY AND CONCLUSIONS 

We present a multi-wavelength analysis of ACT- 
CL J0102-4915 “El Gordo,” a new SZ cluster which 
appears to be the first reported high-redshift Bullet-like clus¬ 
ter undergoing a major merger. We have used a combination 
of optical (VLT), X-ray {Chandra), and infrared ( Spitzer ) 
data to investigate the cluster’s physical properties and find 
that it is an exceptionally massive, X-ray luminous and hot 
system, resulting in the highest SZ signal in the ACT survey. 
From the optical spectra of 89 galaxies confirmed as cluster 
members, we obtain a spectroscopic redshift (z = 0.87008) 
and velocity dispersion (t7 ga i = 1321 ± 106 km s -1 ) for this 
cluster. Our 60-ks Chandra observation on ACIS-I provides 
a spectroscopic temperature of T x = 14.5 ±1.0 keV and 
bolometric luminosity of Lboi = 13.6 x 10 45 h)} } erg s -1 . From 
mass scaling relations for the velocity dispersion. X-ray Y x , 
and SZ distortion, we determine the total mass of the cluster 
to be M 2 oo = (2.16 ±0.32) x 10 15 establishing it as 

the most massive cluster known at z > 0.6. 

From our warm-phase Spitzer IRAC and deep optical imag¬ 
ing, we measure the total stellar mass and constrain the stel¬ 


lar content of the cluster to be <1% of the total mass, in 
broad agreement with other massive clusters. The optical {riz 
bands) and IR (3.6/xm, 4.5/jjn) color-magnitude diagrams for 
spectroscopically-confirmed members show a clearly-defined 
red sequence consistent with passively evolving galaxies, in¬ 
cluding a small fraction of actively star-forming galaxies. 
However the BCG is extremely luminous, significantly bluer 
than the red sequence, and appears as an E+A+[0 II] galaxy 
from the VLT spectra. Such blue BCGs are normally seen 
in the centers of cool core clusters and indeed this BCG sits 
right within the cool, low entropy, enhanced-metal-abundance 
merging core in ACT-CL J0102-4915 revealed by the Chan¬ 
dra data. 

While clusters as massive as ACT-CL J0102-4915 are quite 
rare at its redshift, the cluster does not pose any signifi¬ 
cant difficulty for the standard ACDM cosmology provided 
its mass is in the lower portion of its statistically allowed 
mass range. Our Chandra and VLT observations addition¬ 
ally show that ACT-CL J0102-4915 is undergoing a ma¬ 
jor merger with a mass ratio of approximately 2 to 1 be¬ 
tween its subcomponents. We find no analogous high-mass 
merging systems, with properties broadly similar to ACT- 
CL J0102-4915, within any current large-volume cosmolog¬ 
ical N-body simulations (e.g., MICE, Cubep3m). Thus the 
high mass and merger ratio at this redshift argues for an ex¬ 
ceptional rarity of ACT-CL J0102-4915 within our Universe. 
We expect that more detailed analysis of large cosmological 
simulations and better understanding of the propagation of 
the initial Gaussian fluctuations, particularly at high redshift, 
will be required to compare with the predictions of a ACDM 
model using massive systems like ACT-CL J0102-4915. 

The X-ray surface brightness of our Chandra observation 
reveals a “wake” in the hot gas distribution that we attribute to 
the recent passage of one cluster through the other. The depro- 
jected gas density distribution shows that the wake requires a 
low-density cavity in the interior of the cluster. The timescale 
associated with the refilling of the wake cavity, plus the length 
of the cavity, provides an independent estimate of the speed 
of the merging subcomponent. Under the simple assumption 
that the cavity refills on the sound-crossing time, we estimate 
a time scale of 1 Gyr. The large number of galaxy redshifts 
we measured allows us to identify the two merging compo¬ 
nents and determine their peculiar radial velocity difference 
of around 600 km s _1 . Including this value as the radial com¬ 
ponent of the merger velocity, we estimate both the inclina¬ 
tion angle of the merger axis and the merger speed as 27° and 
1300 km s -1 . An improved estimate of the cavity refilling 
time from N-body hydro simulations should allow for a fairly 
precise estimate of these quantities. The rarity of a merging 
cluster system with the relative masses and velocities of the 
subcomponents of ACT-CL J0102-4915 will require further 
simulations and analysis to determine. 

The next most critical observational step is to obtain strong 
and weak lensing mass estimates and maps from the ground 
and with the Hubble Space Telescope, which we are pursu¬ 
ing. This will allow us to confirm the presence of an offset 
between the baryonic and dark matter distributions, as well 
as provide a robust mass estimate unbiased by merger state 
or cluster astrophysics. A deeper Chandra observation will 
allow better mapping of the spatial temperature distribution 
in the cluster, which may allow identification of shocks in 
the gas. It is now possible to use ALMA to search for these 
same features directly in the SZ effect using the low frequency 
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channel (Band 3 covering 84-116 GHz). A search for a ra¬ 
dio relic is clearly relevant as well since these are typically 
found in merging clusters, although ACT-CL J0102-4915’s 
high redshift may preclude easy detection. As we have out¬ 
lined above, several theoretical steps also need to be pursued 
in order to fully understand the many distinctive features of 
ACT-CL J0102-4915 that we have presented here. We hope 
that this new exceptional cluster will fuel more in-depth sim¬ 
ulations of cluster mergers and motivate continued searches 
for similar high-redshift systems. 
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